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ABSTRACT 

Interferometric  observations  of  two  well-known  Be  stars,  7  Cas  and  (j)  Per,  were  collected  and  analyzed  to 
determine  the  spatial  characteristics  of  their  circumstellar  regions.  The  observations  were  obtained  using  the  Navy 
Prototype  Optical  Interferometer  equipped  with  custom-made  narrowband  filters.  The  filters  isolate  the  Ha  emission 
line  from  the  nearby  continuum  radiation,  which  results  in  an  increased  contrast  between  the  interferometric  signature 
due  to  the  Ha -emitting  circumstellar  region  and  the  central  star.  Because  the  narrowband  filters  do  not  significantly 
attenuate  the  continuum  radiation  at  wavelengths  50  nm  or  more  away  from  the  line,  the  interferometric  signal  in  the 
Ha  channel  is  calibrated  with  respect  to  the  continuum  channels.  The  observations  used  in  this  study  represent  the 
highest  spatial  resolution  measurements  of  the  Ha -emitting  regions  of  Be  stars  obtained  to  date.  These  observations 
allow  us  to  demonstrate  for  the  first  time  that  the  intensity  distribution  in  the  circumstellar  region  of  a  Be  star  cannot 
be  represented  by  uniform  disk  or  ringlike  structures,  whereas  a  Gaussian  intensity  distribution  appears  to  be  fully 
consistent  with  our  observations. 
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1.  INTRODUCTION 

The  application  of  long-baseline  optical  interferometry  to  the 
study  of  classical  Be  stars,  although  still  a  developing  observa¬ 
tional  field,  has  already  resulted  in  significant  contributions  to 
our  understanding  of  these  objects.  This  is  in  part  related  to  the 
fact  that  the  circumstellar  regions  associated  with  the  closest  Be 
stars  can  be  spatially  resolved  using  optical  interferometers  with 
only  modest  baseline  lengths  (10-40  in),  especially  if  the  obser¬ 
vations  are  sensitive  to  the  Ha  line  emission.  The  interferometric 
observations  at  such  baselines  typically  yield  only  angular  sizes 
of  the  emitting  regions  and  in  some  cases  also  the  apparent  ellip- 
ticity.  However,  when  these  types  of  results  are  combined  with 
polarimetry  or  spectroscopy,  a  number  of  different  intrinsic  prop¬ 
erties  of  these  regions  can  then  be  investigated  or  constrained. 

For  example,  Quirrenbach  et  al.  (1997)  combined  optical  in¬ 
terferometric  and  spectropolarimetric  observations  of  seven  Be 
stars  and  showed  that  the  disk  orientations  inferred  from  the  two 
completely  independent  data  sets  agree.  Furthermore,  because 
the  smallest  upper  limit  on  the  disk  opening  angle  derived  by 
Quirrenbach  et  al.  (1997)  from  the  apparent  ellipticity  of  the  cir¬ 
cumstellar  region  was  ~20°,  this  was  the  first  study  that  spatially 
resolved  the  circumstellar  disks  and  supported  the  thin-disk  para¬ 
digm  for  Be  stars.  Wood  et  al.  ( 1 997)  constructed  disk  models  for 
the  Be  star  (  Tau  with  large  (tens  of  degrees)  and  small  (a  few 
degrees)  disk  opening  angles,  in  which  both  models  were  consis¬ 
tent  with  the  spcctropo lari mctiy.  Using  the  interferometric  re¬ 
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suits  of  Quirrenbach  et  al.  (1997),  the  large  opening  angle  solu¬ 
tion  was  rejected  by  Wood  et  al.  (1997). 

Another  example  of  the  synergy  between  spectroscopy  and 
long-baseline  interferometry  was  demonstrated  by  Vakili  et  al. 
(1998),  in  which  the  observations  of  the  Be  star  (  Tau  at  two 
epochs  were  used  to  detect  the  presence  of  a  one-armed  oscilla¬ 
tion  in  the  circumstellar  disk.  In  a  similar  study,  Berio  et  al.  (1999) 
have  shown  that  the  variable  asymmetric  brightness  distribution 
in  the  disk  of  7  Cas  deduced  from  the  optical  interferometric  data 
correlates  with  the  spectral  variations  seen  in  the  H j3  emission 
line,  and  these  in  turn  can  also  be  explained  with  a  precessing  one- 
anned  oscillation  in  the  equatorial  disk.  Although  the  combination 
of  spectroscopic  and  interferometric  observations  is  necessary  in 
studies  related  to  temporal  variability,  these  types  of  data  sets  can 
also  be  combined  to  put  direct  constraints  on  the  physical  condi¬ 
tions  within  the  circumstellar  regions  of  Be  stars.  For  example, 
Tycner  et  al.  (2005)  combined  Ha  emission  profiles  from  spec¬ 
troscopy  with  interferometric  observations  in  Ha  for  a  number 
of  different  Be  stars  and  showed  that  there  exists  a  direct  relation¬ 
ship  between  the  physical  extent  of  the  emitting  region  and  the 
net  Ha  luminosity. 

The  total  number  of  Be  stars  investigated  with  optical  interfer¬ 
ometry  to  date  still  remains  quite  small.  This  can  be  attributed  to 
the  relatively  small  number  of  long-baseline  interferometric  in¬ 
struments  that  are  sensitive  to  one  of  the  strongest  spectral  fea¬ 
tures  produced  within  the  circumstellar  disk,  i.e.,  the  Ha  emission 
line.  Even  the  instruments  that  are  sensitive  to  the  Ha  emission 
line,  or  any  other  emission  line  for  that  matter,  and  have  baselines 
long  enough  to  spatially  resolve  the  circumstellar  regions  still 
need  to  separate  the  signatures  from  the  resolved  disk  and  the  un¬ 
resolved  (or  nearly  unresolved)  central  star  that  is  detected  at 
the  continuum  wavelengths.  For  this  reason,  it  is  imperative  that 
the  instrument  be  configured  such  that  the  contribution  from  the 
stellar  photosphere  to  the  interferometric  observations  at  the  Ha 
line  is  minimized.  This  can  be  accomplished  either  by  sufficiently 
high  spectral  dispersion  or  with  the  use  of  narrowband  filters. 

The  Navy  Prototype  Optical  Interferometer  (NPOI)  is  by  de¬ 
sign  a  multispectral  instrument.  However,  the  continuum  bandpass 
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recorded  by  the  instrument  in  the  Ha  region,  as  defined  by  the 
spectral  channel  characteristics,  can  be  too  wide  for  sources  with 
weak  Ha  emission.  Therefore,  as  part  of  a  proof  of  concept,  we 
equipped  the  NPOl  with  a  set  of  custom  narrowband  filters8  cen¬ 
tered  on  the  Ha  line  to  further  decrease  the  stellar  contribution 
at  the  continuum  wavelengths  in  that  channel  (Pauls  et  al.  2001). 
These  initial  observations  demonstrated  that,  with  only  small 
modifications  to  our  current  instrumental  configuration,  interfer¬ 
ometric  observations  of  Be  stars  can  be  obtained  with  the  narrow- 
band  filters. 

In  this  paper  we  demonstrate  the  quantitative  results  based  on 
the  narrowband  interferometric  observations  of  two  stars  known 
for  their  relatively  strong  Ha  emission,  7  Cassiopeiae  (HR  264) 
and  <f>  Persei  ( HR  496).  Because  the  circumstellar  regions  of  both 
stars  were  resolved  by  long-baseline  interferometry  in  the  past, 
they  are  suitable  targets  for  direct  comparison  of  the  results  pro¬ 
duced  by  the  new  observational  setup  with  those  already  presented 
in  the  literature  (Quirrenbach  et  al.  1997;  Tycner  et  al.  2003).  Also, 
the  interferometric  observations  presented  in  this  study  repre¬ 
sent  some  of  the  highest  spatial  resolution  measurements  of  the 
Ha -emitting  regions  obtained  to  date.  This  allows  us  to  test  the 
various  models  for  intensity  distribution  at  high  enough  spatial 
frequencies  that  the  degeneracy  between  the  various  models  is 
eliminated. 

2.  OBSERVATIONAL  SETUP 

The  NPOl  consists  of  four  stationary  astrometric  stations 
and  a  Y-shaped  imaging  station  configuration  in  which  movable 
siderostat  elements  can  be  repositioned  (see  the  description  of 
the  instrument  by  Armstrong  et  al.  1998).  Currently,  two  mov¬ 
able  stations  are  operational.  The  goal  is  to  have  six  movable  sta¬ 
tions  that  will  allow  simultaneous  observations  with  up  to  six 
reconfigurable  elements.  Although  the  successful  combination 
of  light  from  six  stations  has  already  been  demonstrated  with  the 
NPOl  (Hummel  et  al.  2003),  at  the  time  of  the  observations  pre¬ 
sented  in  this  paper  only  four  stations  were  available  for  obser¬ 
vations.  These  four  stations  gave  access  to  baselines  with  lengths 
in  the  range  of  1 8.9-64.4  m,  as  shown  schematically  in  Figure  1 . 

The  instrument  was  designed  to  accept  stellar  light  from  up 
to  six  siderostats  simultaneously,  in  which  each  element  sends  a 
12  cm  light  beam  through  vacuum  pipes  to  the  beam  combiner 
lab.  At  the  beam  combiner,  the  light  beams  are  split  and  recom¬ 
bined  in  such  a  way  that  the  interference  fringes  between  all  side¬ 
rostat  pahs  (for  a  total  of  15  unique  pairs)  are  constructed.  Figure  2 
illustrates  the  propagation  of  light  for  six  input  beams.  The  three 
output  beams  that  are  intercepted  are  dispersed  by  a  set  of  three 
prisms  onto  individual  lenslet  arrays,  which  in  turn  are  fiber- 
coupled  to  a  cluster  of  photon-counting  avalanche  photodiodes. 
Each  fiber  corresponds  to  a  single  channel  with  spectral  charac¬ 
teristics  set  by  the  position  of  the  lenslet  array  and  the  physical 
width  of  a  single  lenslet  along  the  dispersion  direction. 

The  spectra  of  Be  stars  show  many  spectral  lines  in  emission. 
For  example,  in  the  510-880  nm  region  covered  by  our  interfer¬ 
ometric  observations,  emission  lines  of  elements  such  as  H 1,  He  1, 
O  1,  Si  11,  and  Fe  11  can  be  detected.  However,  except  for  hydro¬ 
gen  lines  all  of  the  emission  lines  are  very  weak,  having  equiv¬ 
alent  widths  (EWs)  of  V0. 1  nm,  and  therefore  are  lost  in  the 
continuum  signal  present  in  each  spectral  channel  (the  widths  of 
the  spectral  channels  range  from  10  to  26  nm).  Even  in  the  case 
of  the  hydrogen  lines  present  in  our  spectral  region  of  interest, 
only  the  Ha  line  has  a  large  enough  EW  to  contribute  a  detect- 
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Fig.  1. —  Schematic  of  the  inner  portion  of  the  NPOl  array.  The  imaging 
(1 circles )  and  astrometric  ( squares )  stations  used  in  the  Ha  observations  pre¬ 
sented  in  this  paper  are  shown  with  filled  symbols.  Baselines  that  could  be 
recorded  on  two  out  of  three  outputs  from  the  beam  combiner  are  also  shown 
with  their  respective  physical  lengths  indicated.  The  baseline  that  was  not  ac¬ 
cessible  to  observations  of  7  Cas  due  to  sky  coverage  limitation  is  shown  with 
the  dashed  line. 

able  interferometric  signal  in  the  Ha  -containing  channel.  There¬ 
fore,  for  the  purpose  of  the  Ha  observations,  the  lenslet  array  is 
aligned  such  that  the  emission  line  is  centered  on  a  single  channel 
that  has  a  spectral  width  of  15  nm.  Because  the  typical  EWs  of 
Ha  emission  lines  in  Be  stars  are  usually  not  more  than  4  nm 
(see,  e.g.,  Table  6  in  Tycner  et  al.  2005),  the  contribution  from  the 
circumstellar  region  to  the  net  signal  is  typically  less  than  25%. 
For  sources  with  very  weak  emission,  the  fractional  contri¬ 
bution  to  the  net  signal  can  be  much  less  than  the  magnitude  of 
the  random  and  systematic  uncertainties,  which  are  typically  at 
the  few  percent  level. 

To  increase  the  contrast  between  the  interferometric  signatures 
from  the  Ha  -emitting  circumstellar  region  and  the  central  star, 
we  have  inserted  a  narrowband  filter  at  each  of  the  three  outputs 
from  the  beam  combiner  (see  Fig.  2).  These  filters  have  a  high 
transmission  (~90%)  in  a  2.8  nm  region  centered  on  the  Ha  line, 
which  suppresses  the  nearby  continuum  emission  but  does  not 
significantly  attenuate  radiation  outside  a  ~  1 00  nm  region  cen¬ 
tered  at  the  line.  Because  we  need  to  detect  the  interferometric 
signal  outside  of  the  spectral  channel  containing  the  Ha  emis¬ 
sion  for  both  calibration  and  fringe  tracking  purposes,  the  filter 
has  high  enough  transmission  in  the  500-600  and  710-850  nm 
regions  that  the  continuum  charnels  are  still  usable.  The  transmis¬ 
sion  curve  of  the  custom-made  filter  is  shown  in  Figure  3  along 
with  the  superposed  locations  of  the  spectral  channels  that  were 
used  to  record  the  fringe  signals. 

In  the  current  instrumental  setup  the  signals  from  up  to  32  chan¬ 
nels  can  be  recorded  simultaneously  by  the  NPOl  electronics. 
Because  these  channels  are  divided  equally  between  two  output 
beams  from  the  beam  combiner,  data  from  one  of  the  three  out¬ 
put  beams  are  currently  not  recorded  (recall  Fig.  2).  This  results 
in  interferometric  observations  that  do  not  sample  all  the  baseline 
configurations  that  are  possible  with  a  given  set  of  elements,  as 
shown  in  Figure  1 .  Furthermore,  because  the  long-delay  lines  are 
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Fig.  2. —  Schematic  of  the  beam  combiner  illustrating  the  propagation  of  light  from  up  to  six  input  beams.  The  three  output  beams  that  are  intercepted  by  the  pick-off 
optics  sample  all  available  baseline  configurations  (up  to  15  element  pairs  with  six  input  beams). 


not  implemented  yet,  the  maximum  path  difference  that  can  be 
introduced  between  a  pair  of  siderostats  (to  compensate  for  the 
different  distances  between  each  station  and  the  star)  is  currently 
35  m,  and  this  results  in  restricted  ranges  of  declination  and  hour 
angle  over  which  a  star  can  be  observed  at  a  given  baseline.  For 
example,  Figure  4  shows  the  accessible  sky  coverage  for  three 
different  configurations;  in  one  case  all  four  elements  (AC,  AW, 
AE,  and  W7)  are  used,  and  in  the  other  two  configurations  ei¬ 
ther  the  AC  or  the  W7  station  is  excluded.  Therefore,  for  7  Cas 
the  observations  utilizing  the  AC-W7  baseline  are  not  currently 
possible,  whereas  for  (j)  Per,  which  is  10°  lower  in  the  sky,  such 
observations  are  possible.  Of  course,  with  the  expected  implemen¬ 
tation  of  the  long-delay  lines,  optical  path  difference  compensa¬ 


tion  of  up  to  '-"-400  m  will  become  possible,  and  the  accessible 
sky  coverage  will  be  significantly  less  restricted. 

3.  OBSERVATIONS 
3.1.  Interferometry 

The  typical  NPOI  observing  and  data  reduction  procedures  are 
described  in  Elummel  et  al.  (2003).  Here  we  provide  a  brief  re¬ 
view  of  the  process.  During  a  30  s  observation  fringe  data  are 
recorded  every  2  ms.  The  squared  visibility  (V2)  estimators  from 
the  2  ms  data  are  averaged  on  1  s  intervals,  and  the  1  s  data  points 
are  processed  using  a  suite  of  custom  reduction  software.  Out¬ 
lier  points  are  flagged  on  the  basis  of  the  residuals  of  the  delay, 
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Fig.  3. — Transmission  curve  for  the  Ho  filter  in  the  spectral  region  covered 
by  the  16  channels  that  were  used  to  record  the  interferometric  signals.  The 
edges  of  the  spectral  channels  are  marked  with  vertical  lines,  and  the  spectral 
regions  for  which  the  signals  were  not  recorded  are  shown  with  hatched  regions. 

seeing  indicators,  photon  rates,  and  squared  visibilities.  Finally, 
depending  on  how  many  1  s  data  points  have  been  flagged  as  bad, 
up  to  30  points  are  averaged  to  produce  a  single  squared  visibility 
measure  for  each  30  s  interval,  also  known  as  a  “scan.”  Typi¬ 
cally,  following  a  procedure  similar  to  that  used  in  obtaining  the 
V 2  data,  the  closure  phase  and  triple  amplitude  can  also  be  ob¬ 
tained  for  each  scan.  However,  because  these  quantities  were  not 
obtained  for  the  data  presented  in  this  paper,  we  do  not  include 
them  in  our  discussion. 

The  V2  value  for  a  single  scan  is  obtained  in  the  same  way  for 
both  the  target  and  the  calibrator  star,  as  well  as  for  the  incoher¬ 
ent  scans  (also  known  as  the  “off-the-fringe  scans”).  The  typical 
observational  sequence  consists  of  a  pair  of  coherent  and  inco¬ 
herent  scans  on  a  target  followed  by  the  same  pair  of  scans  on  a 
calibrator.  The  procedure  is  repeated  for  as  long  as  the  scans  of 
the  target  star  are  acquired.  The  incoherent  scans  are  necessary 
to  estimate  the  additive  bias  temis  affecting  the  squared  visibil¬ 
ity  measures  (see  Hummel  et  al.  2003  for  more  details).  The  bias 
term  is  a  squared  visibility  measure  that  one  obtains  in  the  case  of 
a  completely  incoherent  signal,  which  for  an  ideal  system  should 
be  exactly  zero.  Because  the  calibration  procedure  removes  mul¬ 
tiplicative  effects,  the  subtraction  of  the  additive  bias  temi  must 
be  performed  before  calibration  corrections  are  applied.  We  re¬ 
move  the  bias  terms  from  the  squared  visibilities  of  both  the 
target  and  the  calibrator  stars  following  the  method  described  in 
Wittkowski  et  al.  (2001)  and  Hummel  et  al.  (2003). 

As  we  have  discussed  in  §  2,  7  Cas  could  not  be  observed  on 
the  AC-W7  baseline  because  its  high  declination  places  the  star 
outside  the  accessible  sky  coverage  (defined  by  the  limits  im¬ 
posed  by  the  optical  path  compensation  components).  Having 
slightly  lower  declination,  0  Per  was  not  constrained  by  the  same 
limitation.  Because  the  two  stars  are  separated  in  right  ascension 
by  less  than  50  minutes,  it  was  more  efficient  to  observe  both 
stars  in  the  same  mode  (i.e.,  using  only  three  stations).  The  loss 
of  one  station  for  0  Per  was  compensated  for  by  the  extended 
hour  angle  range  over  which  it  was  possible  to  observe  this  star 
(recall  Fig.  4).  Another  advantage  of  having  only  one  baseline  on 


Hour  Angle  (hr) 

Fig.  4. — Accessible  sky  coverage  for  interferometric  observations  using  dif¬ 
ferent  combinations  of  the  array  elements.  At  a  fixed  declination,  the  largest 
accessible  hour  angle  (HA)  range  is  determined  by  the  limits  imposed  on  the 
siderostats  motion  (to  the  east)  and  on  a  typical  zenith  angle  limit  of  60°  (to  the 
west).  The  accessible  HA  range  is  further  constrained  by  limitations  imposed  by 
optical  path  compensation  components.  The  smallest  area  ( horizontal  hatched 
region)  corresponds  to  observations  at  all  four  stations  (AC,  AE,  AW,  and  W7), 
which  have  a  declination  limit  of  55°.  The  declination  limit,  as  well  as  the  sky 
coverage,  can  be  increased  by  excluding  the  AC  ( diagonal  hatched  region)  or 
the  W7  station  ( vertical  hatched  region).  The  positions  of  the  horizon  ( dashed 
line),  (f>  Per,  and  7  Cas  (dotted  lines)  are  also  indicated. 

each  spectrograph  is  the  elimination  of  multibaseline  cross-talk 
effects  that  may  be  present  in  the  data  (Schmitt  et  al.  2005).  The 
disadvantage  of  having  only  two  simultaneous  baselines  is  that 
some  interferometric  observables,  such  as  closure  phase  and  tri¬ 
ple  amplitude,  could  not  be  obtained  because  they  require  at  least 
three  simultaneous  baselines. 

We  observed  7  Cas  and  </>  Per  during  2004  November  and 
December.  Although  we  lost  many  nights  due  to  poor  weather, 
we  still  managed  to  obtain  ~100  scans  for  each  star  with  almost 
all  the  scans  providing  data  on  two  baselines.  Table  1  shows  a  log 
of  the  observations.  On  2004  December  8  the  vacuum  pipes  be¬ 
tween  the  sixth  and  seventh  station  on  each  of  the  three  arms 
were  disconnected  (see  Fig.  1  for  the  layout  of  the  imaging  array) 
in  order  to  install  safety  vacuum  valves.  This  resulted  in  the  W7 
station  being  unavailable  for  observations  starting  2004  December 
10,  and  therefore,  only  observations  utilizing  the  AE,  AC,  and  AW 
stations  were  recorded  for  the  rest  of  the  observing  run.  On  2004 
December  2  we  temporarily  removed  the  narrowband  filters  from 
the  outputs  of  the  beam  combiner  to  obtain  observations  without 
the  use  of  the  filters  but  with  the  same  instrumental  configuration. 

For  the  purpose  of  the  analysis  presented  in  this  paper,  we  are 
only  interested  in  the  squared  visibilities  from  the  spectral  chan¬ 
nels  that  contain  the  Ha  emission  line.  However,  before  these  V2 
values  can  be  extracted  from  the  observational  data  set,  we  need 
to  calibrate  these  quantities.  We  follow  the  calibration  procedure 
outlined  in  Tycner  et  al.  (2003),  which  we  only  summarize  here. 
The  calibration  procedure  uses  the  continuum  channels  at  a  given 
scan  and  baseline  to  estimate  a  correction  function  that  has  a 
quadratic  dependence  in  A  l,  which  is  then  applied  to  the  squared 
visibility  measure  in  the  spectral  channel  that  contains  the  Ha 
emission  line.  Because  the  correction  function  is  a  slowly  vary¬ 
ing  function  across  channels,  it  cannot  account  for  any  high-order 
chamiel-to-channel  variations.  Therefore,  in  the  current  analysis 
we  estimated  the  V2  amplitudes  of  the  high-order  variations  across 
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TABLE  1 


Observing  Log  for  7  Cas  and  <j>  Per 


UT  Date 

7  Cas 

(Number  of  Scans) 

tj>  Per 

(Number  of  Scans) 

AE-AC 

AE-W7 

AE-AW 

AE-AC 

AE-W7 

AE-AW 

2004  Nov  3 . 

2 

2 

2 

2 

2004  Nov  4 . 

7 

7 

9 

9 

2004  Nov  5 . 

7 

7 

8 

8 

2004  Nov  30 . 

5 

6 

2004  Dec  1 . 

2 

2 

8 

8 

2004  Dec  2a . 

9 

9 

14 

14 

2004  Dec  3 . 

10 

10 

7 

7 

2004  Dec  4 . 

3 

6 

3 

6 

2004  Dec  10 . 

9 

9 

9 

9 

2004  Dec  11 . 

11 

11 

12 

12 

2004  Dec  12 . 

8 

9 

11 

11 

2004  Dec  19 . 

9 

9 

9 

10 

2004  Dec  23 . 

12 

12 

10 

10 

Total . 

49 

40 

98 

51 

51 

112 

Note. — The  AC-W7  baseline  was  not  accessible  to  7  Cas  (see  §  3.1),  and  therefore  the  AC  station  was  not  used  until 
after  2004  December  4,  when  the  W7  station  came  offline  due  to  gate  valve  work. 

“  The  narrowband  Ha  filter  was  not  used  on  2004  December  2. 


the  spectral  charnels  using  the  residuals  that  could  not  be  removed 
with  the  quadratic  polynomials  from  the  scans  of  the  calibrator 
stars.  The  calibrators  used  for  7  Cas  and  0  Per  were  e  Cas  ( HR  542) 
and  C  Cas  (HR  153),  respectively.  Both  calibrators  are  B-type  stars 
and  were  verified  by  spectroscopic  observations  to  have  Ha  in 
absorption.  After  confinning  that  the  channel-to-channel  variations 
established  based  on  observations  of  the  calibrators  are  stable  on 
timescales  longer  than  one  night,  we  calculated  nightly  averages 
that  were  then  divided  out  of  the  scans  of  the  target  stars. 

After  applying  the  bias  and  the  calibration  corrections  to  all 
scans,  we  extracted  the  squared  visibilities  for  the  Ha  channel 
forming  two  large  data  sets,  one  for  7  Cas  (with  169  data  points) 


u  (MX) 

Fig.  5. —  Sampling  of  the  ( u ,  !>)-plane  by  the  Ha  observations  of  7  Cas  on 
three  baselines  with  lengths  of  18.9  (AC-AE),  37.5  (AE-AW),  and  64.4  m 
(AE-W7).  For  comparison,  sample  coverages  (not  limited  by  the  HA  limitations 
shown  in  Fig.  4)  over  6  hr  ranges  in  HA  from  the  meridian  to  the  east  ( dotted 
line)  and  to  the  west  {solid  line)  are  also  shown. 


and  the  other  for  cj)  Per  (with  186  data  points).  The  observations 
obtained  on  2004  December  2,  when  the  narrowband  filters  were 
not  used,  are  treated  as  separate  data  sets.  Figures  5  and  6  show 
the  resulting  («,  t;)-planc  coverage  at  the  AC-AE,  AE-AW,  and 
AE-W7  baselines  in  the  Ha  channel  for  7  Cas  and  (f)  Per,  respec¬ 
tively.  The  corresponding  V 2  data  for  the  two  stars  are  shown  in 
Figures  7  and  8  for  observations  obtained  with  the  narrowband 
filter,  whereas  the  V2  observations  obtained  without  the  filters 
are  shown  in  Figures  9  and  10. 

3.2.  Spectroscopy 

The  interferometric  observations  obtained  with  the  NPOI  do 
not  contain  enough  spectral  information  to  help  us  establish  the 
properties  of  the  Ha  emission  line.  Therefore,  we  observed  7  Cas 
and  4>  Per  using  a  fiber-fed  echelle  spectrograph  (also  known  as 


Fig.  6. —  Same  as  Fig.  5,  but  for  the  (u,  ffl-plane  coverage  of  cj)  Per. 
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Fig.  7. — Calibrated  squared  visibilities  from  the  Ha  channel  of  7  Cas  ob¬ 
tained  at  three  baselines.  The  elliptical  Gaussian  model  ( red  solid  line )  fitted  to 
all  observations  is  shown  at  each  baseline  for  the  same  HA  range  as  defined  by 
the  observations.  The  model  contains  a  contribution  from  the  stellar  photo- 
spheric  disk  ( dashed  line)  that  is  modeled  with  eq.  (2).  Model  curves  evaluated 
at  the  minor  ( top  dotted  line)  and  major  ( bottom  dotted  line)  axes  of  the  elliptical 
Gaussian  model  are  also  shown. 

the  Solar-Stellar  Spectrograph  [SSS]),  which  is  located  at  the 
Lowell  Observatory’s  John  S.  Hall  1.1  m  telescope.  The  SSS  in¬ 
strument  produces  spectra  in  the  Ha  region  with  a  resolving  power 
of  ~  1 0,000  and  a  signal-to-noise  ratio  (S/N )  of  a  few  hundred  in 
the  Ha  region.  The  SSS  observations  were  reduced  using  stan¬ 
dard  reduction  routines  written  in  IDL9  (Hall  et  al.  1 994),  and  the 


9  Interactive  Data  Language  of  RSI,  ITT  Industries,  Boulder,  CO  80301. 


Fig.  9. —  Calibrated  squared  visibilities  from  the  Ha  channel  of  7  Cas  ob¬ 
tained  on  2004  December  2  without  the  use  of  a  narrowband  filter  at  two  base¬ 
lines  (AE-AW  and  AE-W7).  The  best-fit  circularly  symmetric  Gaussian  model 
(i dotted  line)  and  the  stellar  photospheric  disk  ( dashed  line)  are  also  shown. 

resulting  Ha  profiles  of  7  Cas  and  (f)  Per  are  shown  in  Figures  1 1 
and  12,  respectively. 

The  Ha  emission  of  7  Cas  is  known  to  be  stable  on  timescales 
greater  than  1  yr,  and  our  spectra  confirm  this.  The  EWs  of  the 
four  Ha  profiles  shown  in  Figure  1 1  range  from  —  3 .0  to  —  3 .2  nm. 
Because  the  largest  uncertainty  associated  with  these  values  is 
the  precision  of  continuum  normalization,  which  we  estimate  at 
3  %,  for  the  purpose  of  our  study  we  use  an  E  W  of — 3 . 1  ±0.1  nm 
for  the  Ha  emission  of  7  Cas. 


Fig.  10. —  Same  as  Fig.  9,  but  for  the  Ha  observations  of  cf>  Per  with  the  best- 
fit  elliptical  Gaussian  model  shown  ( solid  lines). 


Fig.  8. —  Same  as  Fig.  7,  but  for  the  Ha  observations  of  <j>  Per. 
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Fig.  11. — Ha  profiles  of  7  Cas  obtained  at  four  different  epochs. 


Fig.  12. — Ha  profiles  of  <j)  Per  obtained  at  three  different  epochs. 


As  we  discuss  in  §  5.3,  cj)  Per  is  a  Be+sdO  spectroscopic  binary 
with  a  period  of  127  days  and  is  known  to  show  Ida  variability. 
The  observed  variability  is  attributed  in  part  to  the  radiative  ef¬ 
fects  of  the  hot  secondary  on  the  circumstellar  region  associated 
with  the  primary.  Although  our  temporal  coverage  is  limited,  we 
confirm  the  presence  of  Ela  emission  variability  in  cj)  Per.  The 
EWs  of  the  three  profiles  shown  in  Figure  12  are  —3.8,  —4.2,  and 
—3.6  nm,  which  were  obtained  on  2004  December  3, 2005  March 
3,  and  2005  April  1,  respectively.  Elowever,  because  the  relative 
shape  of  the  Ha  profile  does  not  appear  to  change  significantly, 
we  conclude  that  combining  all  of  the  interferometric  observa¬ 
tions  from  our  observing  run  should  not  result  in  any  significant 
errors  (we  return  to  this  point  in  §  4.2).  For  this  reason  we  adopt 
an  Ha  EW  for  4>  Per  of  —3.9  ±  0.3  nm,  in  which  the  larger 
uncertainty  should  account  for  the  intrinsic  variability  during  our 
interferometric  run. 


4.  THE  ANALYSIS 
4.1.  Models 

The  squared  visibilities  from  the  spectral  channels  containing 
the  Ha  emission  contain  two  signatures,  one  due  to  the  central 
star  and  the  other  due  to  the  circumstellar  region.  Therefore,  these 
data  must  be  modeled  with  a  two-component  model  of  the  form 

^odel=[^^  +  d-^)Fenv]2,  (1) 

where  Vp  and  Vcnv  are  the  visibility  functions  representing  the 
photosphere  of  the  central  star  and  the  circumstellar  envelope, 
respectively,  and  cp  is  the  fractional  contribution  from  the  stellar 
photosphere  to  the  total  flux  in  the  Ha  channel.  Generally,  both 
Vp  and  Vcnv  can  be  complex  functions,  but  because  the  models 
we  consider  in  this  analysis  are  all  point-symmetric,  as  well  as 
concentric,  the  functions  are  real  and  can  be  treated  as  visibility 
amplitudes. 

To  decrease  the  number  of  free  parameters  in  our  models  we 
treat  the  angular  diameter  of  the  central  star  as  a  known  quantity. 
This  is  also  supported  by  the  fact  that  the  central  stars  are  not 
significantly  resolved  even  at  our  longest  baseline.  Furthermore, 


as  we  discussed  in  Tycner  et  al.  (2004),  our  best-fit  models  de¬ 
scribing  the  Ha -emitting  regions  are  not  very  sensitive  to  the 
assumed  diameter  of  the  central  star  (changing  the  diameter  by  a 
factor  of  2  results  in  ~  1%  change  in  the  best-fit  disk  parameters). 
For  the  same  reason  we  can  also  ignore  any  effects  related  to  the 
rapid  rotation  for  which  Be  stars  are  well  known,  because  the  geo¬ 
metric  distortion  and  gravity  darkening  near  the  stellar  equator 
only  weakly  affect  our  best-fit  disk  parameters. 

The  most  widely  used  approach  for  predicting  a  stellar  angular 
diameter  is  to  use  empirically  established  relationships  between 
the  stellar  color  index  and  the  surface  brightness,  which  in  turn  is 
related  to  the  angular  size  (Bames  et  al.  1978;  van  Belle  1999). 
Another  approach  would  be  to  use  tabulations  of  linear  stellar 
diameters  as  a  function  of  spectral  class  (Underhill  et  al.  1979). 
However,  because  there  is  an  intrinsic  scatter  in  stellar  diameters 
at  each  spectral  class  and  the  distance  to  the  source  needs  to  be 
known  to  convert  the  linear  diameter  to  the  angular  size  needed 
for  interferometry,  this  is  not  the  preferred  approach.  Therefore, 
for  7  Cas  and  (f>  Per  we  adopt  the  same  stellar  angular  diameters 
that  were  used  by  Quirrenbach  et  al.  (1997),  which  were  derived 
using  the  photometric  relations  of  Bames  et  al.  (1978).  The  angu¬ 
lar  diameters  we  adopt  for  7  Cas  and  0  Per  are  0.56  and  0.39  mas, 
respectively. 

We  model  the  stellar  photosphere  of  each  star  with  a  circularly 
symmetric  unifomi  disk  brightness  distribution.  This  ensures  that 
the  central  star  is  modeled  in  exactly  the  same  way,  at  the  contin¬ 
uum  charnels  during  the  calibration  procedure  (see  Tycner  et  al. 
2003  for  more  details)  and  at  the  Ha  channel  via  equation  (1). 
Therefore,  we  represent  the  visibility  function  of  the  photospheric 
component  of  an  angular  diameter  a  with 


VJu,  v)  =  2 


J\  (7 xa\J  u1  +  v2) 
itas/ur  +  v2 


(2) 


where  J\  is  a  first-order  Bessel  function  and  u  and  v  are  the  spatial 
frequencies,  which  are  given  by  the  east-west  and  south-north 
components  of  the  projected  baseline  on  the  plane  of  the  sky 
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TABLE  2 

Best-Fit  Model  Parameters  for  7  Cas  and  <f>  Per 


6  0 

Model  (mas)  t  r  (deg)  cp  xl 


7  Cas 

UD .  4.72  ±  0.03  ...  0.89  ±  0.02  29.6  ±  2.5  0.53  ±  0.01  6.0 

UD* .  6.49  ±  0.06  ...  0.52  ±  0.03  34.7  ±  1.5  0.59  ±  0.01  1.6 

UR .  7.91  ±  0.09  0.04-0.90  0.58  ±  0.03  5-40  0.58  ±0.13  5.7 

UR* .  5.60  ±  0.53  0.77  ±  0.16  0.52  ±  0.03  36.0  ±  1.5  0.67  ±  0.01  1.7 

GD .  3.59  ±  0.04  ...  0.58  ±  0.03  31.2  ±1.2  0.51  ±  0.01  1.4 

GD* .  3.64  ±  0.05  ...  0.56  ±  0.03  30.6  ±1.4  0.51  ±  0.01  1.4 


cf)  Per 

UD .  5.43  ±0.13  ...  0.27  ±0.01  -61.5  ±  0.6  0.52  ±  0.01  1.4 

UR .  5.23  ±  0.33  0.58  ±  0.12  0.26  ±  0.01  -61.7  ±  0.6  0.62  ±  0.01  1.4 

GD .  2.89  ±  0.09  ...  0.27  ±  0.01  -61.5  ±  0.6  0.39  ±  0.02  1.4 


Notes. — For  the  UD  and  UR  models  9  corresponds  to  the  angular  diameter  of  the  major  axis,  and  for  the  GD  model  it  corresponds 
to  the  FWHM  of  the  major  axis  of  the  elliptical  Gaussian.  Asterisks  denote  models  fitted  to  data  from  the  shortest  two  baselines  only. 


divided  by  wavelength  (see  §  4.1  in  Thompson  et  al.  [2001]  for 
the  definition  of  the  [«,  i>]-plane). 

Our  interferometric  observations  were  obtained  at  high  spatial 
frequencies  (i.e.,  at  high  spatial  resolution),  and  therefore,  we 
can  compare  different  models  of  brightness  distribution.  Previ¬ 
ously  published  work  on  disk  models  of  Be  stars  has  demon¬ 
strated  that  the  thermal  structure  of  the  circumstellar  disks  can  be 
quite  complex  (Millar  &  Marlborough  1998;  Millar  et  al.  2000; 
Carciofi  &  Bjorkman  2006).  However,  for  this  first  investigation 
with  this  new  observational  technique  we  choose  to  model  our 
data  with  three  simple  models:  a  uniform  disk  (UD),  a  uniform 
ring  (UR),  and  a  Gaussian  distribution  (GD).  The  visibility  am¬ 
plitudes  for  all  three  models  can  be  written  as 


Va 


2Ji(7r0UDs) 

tt9vds 

2(1  -  e2)-1 


Ji^urs) 


exp 


(tt9gds)2 


4  In  2 


e2Ji(ire9VRs) 

TreduRS 


UD, 

UR,  (3) 
GD, 


where  9\ jD  and  #UR  correspond  to  the  angular  diameters  of  the 
UD  and  UR  models,  respectively,  and  (9GD  corresponds  to  the  full 
width  at  half-maximum  (FWHM)  of  the  Gaussian  model.  Be¬ 
cause  we  allow  all  three  models  to  have  elliptical  distribution  on 
the  sky,  all  of  the  above  diameters  describe  the  dimensions  along 
the  major  axis.  The  other  variables  in  equation  (3)  are  e,  which 
describes  the  inner  diameter  of  the  ring  model  along  the  major 
axis  (in  units  of  d(;R ),  and  s,  which  is  given  by 


s  =  \J r2(u  cos  4>  —  v  sin  4>)2  +  (u  sin  0  +  v  cos  0)2, 


(4) 


where  r  is  the  axial  ratio  and  <f>  is  the  position  angle  (P.A.;  mea¬ 
sured  east  from  north)  of  the  major  axis  (when  0  <  r  <  1).  For 
circularly  symmetric  structures  (r  =  1 ),  equation  (4)  reduces  to  a 
simple  expression  for  a  radial  spatial  frequency  in  the  («,  tj-planc 
given  by  (u2  +  v2)V2. 


4.2.  Model  Fitting 

To  obtain  a  best  fit  for  each  of  the  three  models,  we  used  a 
nonlinear  least-squares  IDL  procedure  based  on  the  Levenberg- 


Marquardt  method  (Press  et  al.  1992).  Each  model  is  represented 
by  equation  ( 1 ),  with  the  appropriate  expression  for  the  envelope 
component  from  equation  (3),  and  a  fixed  visibility  function  for 
the  photospheric  component  (eq.  [2]  with  a  fixed  stellar  diameter 
a).  Therefore,  the  UD  and  GD  models  have  four  free  parameters, 
cp,  9,  r,  and  q !>,  whereas  the  UR  model  has  e  as  an  extra  parameter. 
Table  2  lists  the  best-fit  values  for  the  parameters  for  all  three 
models  along  with  the  reduced  x2  values.  By  inspecting  the  re¬ 
duced  x2  values  in  Table  2  for  the  three  7  Cas  models,  it  is  evi¬ 
dent  that  the  UD  and  UR  models  cannot  represent  the  observations 
obtained  at  all  three  baselines.  Therefore,  we  have  also  fitted  mod¬ 
els  to  observations  of  7  Cas  from  the  shortest  two  baselines  only, 
at  which  all  three  models  yield  similar  quality  fits  (these  models 
are  also  listed  in  Table  2).  To  check  for  any  possible  variations  in 
the  models  that  could  be  attributed  to  the  temporal  variability  of 
our  sources  (especially  in  the  case  of  c j>  Per),  we  fitted  models  to 
various  subsets  of  our  large  data  set.  Although  typically  the  un¬ 
certainties  of  the  best-fit  parameters  are  larger  for  models  fitted  to 
smaller  data  sets,  all  results  were  consistent  with  our  solutions 
based  on  all  observations. 

To  illustrate  why  the  UD  and  UR  models  fail  to  represent  the 
observations  of  7  Cas  obtained  at  all  three  baselines,  in  Figure  1 3 
we  plot  the  best-fit  UD,  UR,  and  GD  models,  which  were  only 
fitted  to  data  from  the  shortest  two  baselines.  Each  baseline  pro¬ 
vides  data  over  a  different  spatial  frequency  range  (because  of 
different  physical  lengths),  and  these  ranges  are  marked  with  thick 
solid  lines  in  Figure  13.  A  small  fraction  of  each  range  corre¬ 
sponds  to  an  instant  at  which  the  projected  baseline  on  the  sky 
was  oriented  along  the  major  axis  for  which  these  model  curves 
were  evaluated.  The  model  curves  demonstrate  the  inherent  de¬ 
generacy  present  between  the  three  models  at  low  spatial  fre¬ 
quencies  (i.e.,  at  low  spatial  resolutions  the  three  models  look 
the  same).  However,  the  longest  baseline  provides  data  at  high 
enough  spatial  frequencies  that  the  degeneracy  between  the  three 
models  is  eliminated  (see  also  Fig.  7).  We  conclude  that  the  UD 
and  UR  models  are  inconsistent  with  the  data  obtained  at  the 
longest  baseline  (largest  spatial  frequencies),  whereas  the  GD 
models  fit  the  observations  at  all  baselines. 

The  observations  of  4>  Per  are  more  challenging  because  the 
angular  size  of  its  circumstellar  region  is  smaller  than  for  7  Cas, 
and  the  projected  baselines  on  the  plane  of  the  sky  are  not  ori¬ 
ented  along  its  major  axis.  This  results  in  interferometric  obser¬ 
vations  that  do  not  resolve  the  region  sufficiently  well  to  break 
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Fig.  13. — Model  squared  visibilities  for  the  uniform  disk  ( solid  line),  ring 
(i dash-dotted  line),  and  Gaussian  ( dotted  line )  models  of  7  Cas.  All  three  models 
contain  a  contribution  from  a  central  star  ( dashed  line),  which  is  modeled  using 
eq.  (2),  and  are  evaluated  along  their  major  axes  (i.e.,  along  their  largest  angular 
extent).  The  spatial  frequency  ranges  sampled  by  the  three  baselines  are  indi¬ 
cated  with  thick  solid  lines.  The  models  shown  were  fitted  to  the  data  at  the  two 
shortest  baselines  only. 


Fig.  14. —  SameasFig.  13,  but  for  models  based  on  observations  of  <j>  Per  and 
evaluated  at  three  different  orientations  with  respect  to  the  major  axis.  The  three 
families  of  curves  correspond  to  1 00%  ( bottom  three  curves),  50%  ( middle  three 
overlapping  curves),  and  27%  ( top  three  overlapping  curves)  of  the  size  of  the 
major  axis.  Because  the  observations  do  not  sample  the  disk  along  the  major  axis, 
the  spatial  frequency  ranges  sampled  by  the  three  baselines  are  indicated  only  along 
the  minor  axis  and  where  the  extent  is  50%  of  the  major  axis  {thick  solid  lines). 


the  degeneracy  between  the  different  models.  This  is  also  the 
reason  that  all  three  models  for  1 j)  Per  listed  in  Table  2  yield  the 
same  xl  values.  Figure  14  demonstrates  the  differences  between 
the  model  curves  based  on  the  three  best-fit  models  evaluated  at 
three  different  projections.  The  three  projections  were  chosen  to 
correspond  to  the  minor  axis  (that  is,  27%  of  the  major  axis),  to  a 
direction  with  an  extent  of  50%  of  the  major  axis,  and  to  the  ma¬ 
jor  axis.  Because  we  do  not  have  observations  that  resolve  the 
major  axis  (recall  Fig.  8),  the  spatial  frequency  ranges  over  which 
we  have  data  are  marked  only  on  the  curves  that  correspond  to  the 
minor  axis  and  50%  of  the  major  (Fig.  14,  thick  solid  lines).  Based 
on  Figure  14  we  conclude  that  in  order  to  break  the  degeneracy 
between  the  models,  observations  either  at  higher  spatial  frequen¬ 
cies  or  along  the  major  axis  are  required  (as  demonstrated  in  the 
case  of  7  Cas). 


4.3.  Gaussian  Distribution 


For  0  Per  we  cannot  exclude  the  UD  and  UR  models  in  favor 
of  the  GD  model  based  on  interferometric  data  alone;  however, 
we  can  demonstrate  that  the  GD  model  is  the  preferred  solution. 
This  is  because  the  three  models  that  fit  the  4>  Per  dafa  equally 
well  have  different  best-fit  values  for  the  cp  parameter,  which  can 
be  constrained  by  spectroscopy.  For  example,  if  the  EW  of  the 
continuum  flux  from  the  star  in  the  Fla  channel  is  F,  (which  has 
units  of  length)  and  the  EW  of  the  net  Ha  emission  is  ^Ha,  then 
the  fractional  photospheric  contribution  in  the  absence  of  a  nar¬ 
rowband  filter  (c*)  can  be  expressed  as 


F* 


Cp  F*  + 


Ha 


(5) 


With  the  narrowband  filter,  F,  gets  reduced  to  tF, ,  where  t  is  the 
fractional  transmission  in  the  Ha  channel.  The  Ha  emission  gets 
reduced  by  the  filter  as  well,  but  because  the  Ha  emission  lines 


have  FWHM  A 1  nm,  which  is  less  than  the  width  of  2.8  nm  of 
the  narrowband  region,  we  can  approximate  the  transmission  at 
the  Ha  line  (//)  by  the  peak  transmission  of  92%  of  the  narrow- 
band  region  (recall  Fig.  3).  Therefore,  the  cp  parameter  for  the 
case  including  the  filter  can  be  written  as 

_  tF ; 

Cp  tF *  +  VEHa  '  1  j 

Taking  the  ratio  of  equations  (5)  and  (6)  allows  us  to  obtain  an 
expression  for  t  of  the  form 


t  -  •  // 


!/<•;  - 1 

1  /cp  - 1  ’ 


(?) 


which  does  not  depend  on  7%,  or  ^Ha- 

Equation  (7)  requires  the  photospheric  contributions  cp  and  c* 
to  be  known.  For  cp  we  use  the  values  listed  in  Table  2  obtained 
for  the  GD  model  fits  to  observations  with  the  narrowband  filters. 
To  obtain  the  corresponding  values  for  c*,  we  fit  GD  models  to 
observations  obtained  without  the  filters  (i.e.,  the  observations 
obtained  on  2004  December  2).  The  squared  visibilities  for  the 
best-fit  models  to  the  2004  December  2  data  are  shown  in  Figures  9 
and  10,  for  7  Cas  and  0  Per,  respectively.  In  both  cases  the  best-fit 
GD  model  parameters  (except  for  cp)  agree  with  the  parame¬ 
ters  listed  in  Table  2,  even  for  7  Cas,  for  which  we  did  not  have 
enough  points  to  constrain  an  elliptical  model.  These  fits  yield  c* 
values  of  0.85  and  0.78  for  7  Cas  and  (j)  Per,  respectively.  Using 
equation  (7)  with  the  appropriate  values  for  7  Cas  and  (j)  Per  we 
obtain  the  same  estimate  for  t  of  0.1 7  ±  0.01.  It  is  interesting  to 
note  that  this  1 7%  transmission  is  fully  consistent  with  the  de¬ 
tected  drop  in  photon  counts  in  the  Ha  channel  when  the  nar¬ 
rowband  filters  are  introduced  for  observations  of  calibrator  stars 
(which  do  not  have  Ha  in  emission). 
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TABLE  3 


ClRCUMSTELLAR  REGIONS  OF  7  CaS  AND  <j)  Per 


Description 

Symbol 

7  Cas 

< j>  Per 

Disk  size  (mas) . 

$GD 

3.59  ±  0.04 

2.89  ±  0.09 

Axial  ratio . 

r 

0.58  ±  0.03 

0.27  ±  0.01 

RA.  of  major  axis  (deg) . 

<t> 

31.2  ±  1.2 

-61.5  ±  0.6 

Inclination  (deg) . 

i 

>55 

>74 

Half-opening  angle . 

<17 

<8 

Because  our  estimate  for  t  does  not  depend  on  F,  or  EUa,  we 
can  use  equation  (6)  to  obtain  a  value  for  cp  based  on  spectroscopic 
observations.  We  approximate  F,  with  the  spectral  width  of  the 
Ha  channel  (which  has  been  measured  to  be  15  ±  1  nrn).  To  ob¬ 
tain  the  EW  of  the  net  Ha  emission  for  4>  Per,  we  use  the  Ha  EW 
from  §  3.2  and  add  a  small  contribution  (0.37  nm)  due  to  the  ab¬ 
sorption  line  that  has  been  filled  in  (using  the  same  procedure  as 
outlined  in  §  4.2  of  Tycner  et  al.  2005).  The  net  Ha  emission  of 
4.27  ±  0.30  nm  for  0  Per  results  in  cp  of  0.39  ±  0.03.  This  es¬ 
timate  is  significantly  lower  than  the  values  listed  in  Table  2  for 
the  UD  and  UR  models,  which  have  best-fit  cp  values  of  0.52  and 
0.62,  respectively,  but  it  is  fully  consistent  with  the  0.39  ±  0.02 
value  obtained  for  the  GD  model.  We  interpret  this  as  an  indi¬ 
cation  that  the  brightness  distribution  of  the  circumstellar  disk  of 
4>  Per  is  best  represented  by  a  GD  model.  It  is  interesting  to  note 
that  if  we  follow  the  same  reasoning  for  7  Cas  data  from  the  short¬ 
est  two  baselines  (for  which  the  degeneracy  between  the  three 
models  also  exists),  we  arrive  at  the  same  conclusion  that  the 
GD  model  is  the  preferred  solution.  For  7  Cas,  based  on  its  Ha 
EW  of  —3.1  ±  0.1  nm  and  an  estimated  component  account¬ 
ing  for  the  absorption  line  of  0.26  nm,  we  estimate  that  cp  = 
0.45  zb  0.03,  which  again  is  closest  to  the  value  obtained  for  the 
GD  model. 

We  conclude  that  out  of  the  three  models  we  have  considered 
in  our  analysis,  the  GD  model  is  the  only  model  that  fits  7  Cas 
observations  and  is  the  preferred  model  for  0  Per  based  on  spec¬ 
troscopic  constraints.  Therefore,  in  Table  3  we  list  our  final  best- 
fit  parameters  describing  the  circumstellar  regions  of  7  Cas  and 
<f>  Per  based  on  the  GD  model  fits.  The  model  squared  visibilities 
calculated  based  on  these  models  are  shown  with  red  solid  lines 
in  Figures  7  and  8  for  7  Cas  and  0  Per,  respectively. 

5.  DISCUSSION 
5.1.  Circumstellar  Region 

Because  the  number  of  published  parameters  describing  the 
circumstellar  regions  of  7  Cas  and  (j>  Per  is  still  very  small,  it  is 
worthwhile  to  compare  our  results  to  those  published  in  the  lit¬ 
erature.  Both  Quirrenbach  et  al.  (1997)  and  Tycner  et  al.  (2003) 
fitted  a  GD  model  to  their  observations,  and  therefore,  their  model 
values  can  be  compared  directly  with  those  listed  in  Table  3 . 

Using  observational  data  from  the  Mark  III  interferometer, 
Quirrenbach  et  al.  ( 1 997)  obtained  best-fit  model  parameters  for 
7  Cas  of  8  =  3.47  zb  0.02  mas,  r  =  0.70  zb  0.02,  and  P.A.  = 
19°  zb  2°.  Tycner  et  al.  (2003)  reported  results  based  on  older 
NPOI  observations  (with  a  maximum  baseline  of  37.5  m),  which 
had  8  =  3.67  zb  0.09  mas,  r  =  0.79  zb  0.03,  and  P.A.  =  32°  zb 
5° .  Our  new  determination  of  the  angular  size  of  the  major  axis  as 
3.59  zb  0.04  fully  confirms  the  previous  estimates.  The  apparent 
axial  ratio  of  0.58  zb  0.03  that  we  find  for  7  Cas  is  smaller  than 
the  published  values;  this  might  be  related  to  the  fact  that  our 
interferometric  observations  cover  a  much  larger  range  of  base¬ 
line  projection  angles  on  the  plane  of  the  sky  (recall  Fig.  5).  Our 


best-fit  value  for  the  P.  A.  of  3 1  °2  zb  1  °2  agrees  with  our  previous 
estimate  (Tycner  et  al.  2003)  but  differs  from  the  P.A.  reported  by 
Quirrenbach  et  al.  ( 1 997)  whose  value  was  at  a  right  angle  to  the 
polarization  vector  that  they  obtained  from  polarimetry.  Possible 
sources  for  this  discrepancy  could  be  the  intrinsic  variability  of 
the  polarization  of  the  source  (see,  e.g.,  Table  6  in  Quirrenbach 
et  al.  1997),  the  residual  effects  associated  with  the  removal  of 
the  interstellar  polarization,  or  the  effects  of  a  nonaxisymmetric 
scattering  surface  (i.e.,  only  for  axisymmetric  sources  can  one 
expect  the  polarization  vector  to  be  exactly  perpendicular  to  the 
plane  of  the  disk). 

For  4>  Per,  Quirrenbach  et  al.  (1997)  obtained  model  values  of 
8  =  2.67  ±  0.20  mas,  r  =  0.46  ±  0.04,  and  P.A.  =  -62°  ±  5°, 
but  they  point  out  that  the  value  for  r  might  be  an  upper  limit  be¬ 
cause  their  baselines  were  not  long  enough  to  resolve  the  minor 
axis.  Indeed,  our  results  based  on  observations  that  use  baselines 
that  are  more  than  twice  as  long  as  those  used  in  their  study  yield 
a  best-fit  value  of  0.27  ±  0.01  for  the  axial  ratio.  However,  the 
angular  size  of  the  major  axis  of  2.89  ±  0.09  mas  and  P.A.  of 
— 6 1  ?5  ±  0)6  that  we  obtain  in  this  study  fully  agree  with  their 
values.  Since  (j>  Per  was  not  observed  by  any  other  interferometer 
than  Mark  III,  our  results  represent  the  first  confirmation  of  the 
values  reported  by  Quirrenbach  et  al.  (1997). 

5.2.  Inclination  and  Disk  Opening  Angles 

Our  interferometric  observations  of  the  spatially  resolved  cir¬ 
cumstellar  regions,  which  show  an  apparent  ellipticity,  can  be  used 
to  estimate  the  inclination10  and  disk  opening  angles.  For  example, 
if  we  assume  that  the  circumstellar  disk  is  circularly  symmetric, 
we  can  obtain  a  lower  limit  on  the  inclination  angle  (i)  with  the 
observed  axial  ratio  (r),  using 

i  ±  cos-1  r.  (8) 

The  minimum  value  corresponds  to  a  geometrically  thin  disk,  in 
which  the  entire  signature  of  the  apparent  ellipticity  is  interpreted 
as  a  projection  effect  in  this  case.  Similarly,  if  we  assume  that  the 
geometry  of  the  circumstellar  region  can  be  represented  with  a 
simple  equatorial  disk  model  with  a  half-opening  angle  Hg  and 
radius  7?disk  (sec,  e.g.,  Fig.  2  in  Waters  et  al.  1987),  we  can  then 
obtain  an  upper  limit  on  Hg  using 

Hg<  sin  \r/2\  (9) 

where  the  maximum  value  corresponds  to  a  system  viewed 
edge-on. 

In  Table  3  we  list  the  estimated  limits  on  the  inclination  and 
half-opening  angles  for  both  7  Cas  and  0  Per  based  on  their  best- 
fit  values  of  r  obtained  from  the  elliptical  Gaussian  models.  In¬ 
terestingly,  for  (j)  Per  the  lower  limit  of  ~74  for  i  is  consistent 
with  the  80° -88°  range  originally  proposed  by  Poeckert  (1981). 
Although  Poeckert  (1981)  did  not  describe  the  circumstellar 
region  in  terms  of  an  opening  angle,  the  disk  thickness  of  (1.3- 
2.7)7?*  and  the  disk  radius  of  7.77?*  can  be  translated  to  a  half¬ 
opening  angle  between  5°  and  10°,  which  again  is  consistent  with 
our  upper  limit  of  ^8°. 

It  is  instructive  to  compare  the  upper  limits  we  obtain  for  Hg 
of  7  Cas  and  4>  Per  to  the  value  of  15°  adopted  by  Waters  et  al. 
(1987)  in  their  study  of  almost  60  Be  stars  based  on  far-IR  IRAS 
observations.  Waters  et  al.  (1987)  estimated  a  number  of  differ¬ 
ent  Be  stellar  characteristics,  such  as  the  mass-loss  rate  and  the 
disk  density,  that  had  a  dependence  on  Hg .  Although  the  authors 


10  The  angle  between  the  normal  to  the  plane  of  the  disk  and  the  line  of  sight. 
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assumed  an  error  of  a  factor  of  1 .5  in  sin  Hg ,  our  upper  limit  for 
the  disk  half-opening  angle  for  0  Per  is  half  their  adopted  value. 
In  the  case  of  0  Per,  the  mass-loss  rate  derived  by  Waters  et  al. 
(1987)  is  overestimated  by  at  least  a  factor  of  2,  or  more  if  the 
source  is  not  viewed  edge-on.  Our  upper  limit  on  the  half-opening 
angle  of  17°  for  7  Cas  is  consistent  with  the  value  adopted  by 
Waters  et  al.  (1987).  However,  based  on  the  line  profile  classi¬ 
fication  scheme  established  by  Hanuschik  et  al.  (1996)  we  can 
conclude  that  7  Cas  is  not  viewed  edge-on,  and  therefore  its  half¬ 
opening  angle  is  most  likely  less  than  17°.  In  the  case  of  the  Be 
star  (  Tau,  for  which  the  apparent  ellipticity  of  the  circumstellar 
region  has  also  been  detected  (Quirrenbach  et  al.  1997;  Tycner 
et  al.  2004),  one  can  conclude  that  Hg^  9° .  It  is  possible  that 
other  Be  stars  also  have  small  disk  opening  angles.  To  further  test 
the  generality  of  this  hypothesis,  future  interferometric  observa¬ 
tions  should  concentrate  on  systems  that  are  thought  to  be  viewed 
at  large  inclination  angles  (e.g.,  using  the  Ha  profile  classification 
of  Hanuschik  et  al.  1 996). 

We  can  also  compare  our  interferometrically  determined  open¬ 
ing  angles  with  those  predicted  by  Stee  (2003),  who  has  calculated 
the  disk  opening  angles  for  a  number  of  Be  stars  based  on  the  flux 
in  the  Brackett  continuum  near  2.2  /n m  using  the  SIMECA  code 
(Stee  &  de  Araujo  1 994;  Stee  et  al.  1 995).  The  main  characteristic 
of  the  SIMECA  model  is  that  for  a  fixed  disk  density  at  the  stellar 
photosphere,  the  opening  angle  is  directly  proportional  to  the 
continuum  flux,  so  that  a  large  flux  in  the  K  band  corresponds  to  a 
large  opening  angle.  For  example,  Stee  (2003)  obtained  a  value 
of  1 4°  for  the  full  opening  angle  of  the  disk  of  0  Per,  which  agrees 
with  our  upper  value  of  1 6° .  However,  our  upper  limit  for  the  full 
opening  angle  for  7  Cas  is  34° ,  which  is  significantly  less  than  the 
54°  obtained  by  Stee  (2003)  for  this  star.  This  discrepancy  may 
be  due  to  the  fact  that  Stee  (2003)  might  have  used  ap¬ 
parent  and  not  absolute  magnitudes  to  derive  the  disk  opening 
angles,  in  which  case  the  values  reported  in  that  study  should  be 
reevaluated. 

5.3.  Disk  Truncation 

Tycner  et  al.  (2004)  showed  that  the  Ha  -emitting  region  of  the 
binary  Be  star  C,  Tau  was  well  within  the  estimated  Roche  lobe  of 
the  primary,  which  suggests  that  the  disk  is  tmncated.  Another 
example  of  a  possible  disk  truncation  was  presented  by  Chesneau 
et  al.  (2005)  for  a  Ara,  for  which  the  Very  Large  Telescope  Inter¬ 
ferometer  operating  in  the  N  band  did  not  resolve  any  structure 
and  therefore  allowed  the  authors  to  put  an  upper  limit  on  the 
extent  of  the  circumstellar  region.  Because  this  limit  was  smaller 
than  the  expected  spatial  extent  based  on  the  models  that  were 
fitted  to  the  Balmer  emission  lines,  the  authors  suggest  that  disk 
truncation  by  an  unseen  companion  might  be  occurring. 

Evidence  is  accumulating  that  suggests  that  7  Cas  is  a  member 
of  a  binary  system.  Hannanec  et  al.  (2000)  published  the  first 
orbital  radial  velocity  curve  for  7  Cas  based  on  observations 
collected  in  the  628-672  nm  range  from  1 993  to  2000.  Their  anal¬ 
ysis  suggests  that  this  star  is  the  primary  component  of  a  spec¬ 
troscopic  binaiy  with  a  period  of  ^204  days  and  an  eccentricity  of 
0.26.  This  result  was  confirmed  by  Miroshnichenko  et  al.  (2002) 
using  high-resolution  spectroscopic  observations  of  the  Ha  line 
obtained  over  a  similar  time  period.  They  report  a  periodic  change 
of  205  days  in  this  line,  which  they  also  attribute  to  the  binary 
system.  Although  Hannanec  et  al.  (2000)  conclude  that  the  sec¬ 
ondary  could  be  either  a  hot  compact  object  or  a  low-luminosity 
late-type  star,  they  estimate  that  the  separation  between  the  com¬ 
ponents  at  periastron  could  be  between  250  and  300  RQ.  The 
distance  to  7  Cas  is  188  ±  20  pc  based  on  Hipparcos  satellite 
measurements  (Perryman  et  al.  1997),  which  places  the  separa¬ 


tion  between  the  components  at  6.2-7.4  mas.  Our  best-fit  value 
for  the  FWHM  of  the  major  axis  is  3.59  ±  0.04,  and  there¬ 
fore,  the  circumstellar  disk  is  well  within  the  orbit  established  by 
Hamianec  et  al.  (2000).  However,  because  the  binary  parameters 
of  7  Cas  are  not  well  established  this  comparison  is  tentative. 

The  binary  nature  of  0  Per  is  much  better  established,  since 
it  was  recognized  as  a  binary  in  the  early  1900s  (Cannon  1910; 
Ludendorff  1911).  Poeckert  (1981)  documents  the  development 
in  the  understanding  of  c f i  Per ’s  binary  nature  over  the  succeeding 
decades,  and  suggests  that  a  subdwarf  O  star  is  the  secondary. 
Gies  et  al.  (1998)  detected  the  secondary  in  ultraviolet  spectra  ob¬ 
tained  with  the  Hubble  Space  Telescope  and  confirmed  Poeckert’s 
prediction.  Along  with  HR  2142  (Waters  et  al.  1991),  59  Cyg 
(Maintz  2003),  and  FY  CMa  (Rivinius  et  al.  2004),  0  Per  is  now 
part  of  a  group  of  four  candidate  Be+sdO  binaries.  Be+sdO  bi¬ 
naries  may  result  from  a  spin-up  of  the  B  star  as  a  result  of  mass 
transfer  from  the  progenitor  of  the  evolved  subdwarf  companion. 

Estimates  of  the  orbital  parameters  of  0  Per  are  available  in 
Bozic  et  al.  (1995)  and  Gies  et  al.  (1998).  The  semimajor  axes 
of  these  two  orbital  solutions  range  from  approximately  230  to 
310  R..  Using  the  Hipparcos  distance  of  220+^  pc  for  0  Per 
(Perryman  et  al.  1997),  we  expect  an  angular  separation  of  4.9- 
6.6  mas  between  the  components  of  0  Per.  As  we  would  expect, 
this  separation  is  larger  than  the  angular  radius  of  any  of  the  mod¬ 
els  of  0  Per  listed  in  Table  2.  The  radial  extent  of  the  Gaussian 
model  can  be  approximated  with  the  FWHM  measure  of  3.12  ± 
0.08  mas  (148  ±  30  R.)  we  obtain  for  0  Per  (i.e.,  assuming  that 
the  radial  extent  is  twice  the  half-width  of  the  Gaussian).  This 
disk  radius  is  very  close  to  the  178-204  RQ  Roche  radius  of  the 
primary  obtained  by  Bozic  et  al.  (1995).  This  is  most  likely  an¬ 
other  example  of  a  Ha  -emitting  region  that  is  close  in  size  to  the 
actual  extent  of  the  circumstellar  region,  and  this  suggests  that 
disk  truncation  is  occurring  in  this  system. 

We  should  note  that  in  the  case  of  0  Per  the  presence  of  a 
truncated  disk  has  been  predicted.  Waters  (1986)  used  the  IRAS 
observations  of  0  Per  at  12, 25,  and  60  pi n  to  measure  the  IR  ex¬ 
cess  and  constrained  the  density  distribution  of  p(R)  =  p0(r/R*)~" 
with  n  =  3.1  for  the  disk  model.  Because  he  obtained  a  value  of 
n  —  2.4  for  two  other  Be  stars,  6  Cen  and  y  Oph,  which  also  hap¬ 
pens  to  correspond  to  a  velocity  law  that  has  a  more  gradual  in¬ 
crease  with  distance  (something  that  is  predicted  based  on  the  Ha 
spectroscopy),  he  argued  that  a  density  distribution  with  n  =  2.4 
might  be  more  appropriate  for  0  Per,  which  could  be  only  sat¬ 
isfied  if  the  disk  model  was  truncated  at  ~6.5R*  (or  ~46  R  -). 
The  apparent  discrepancy  between  our  values  for  the  size  of  the 
truncated  disk  might  be  related  to  the  different  wavelength  re¬ 
gimes  of  our  studies,  as  well  as  to  the  assumption  made  by  Waters 
(1986)  that  the  system  is  viewed  pole-on,  whereas  i  might  be  very 
close  to  90°. 

Poeckert  (1981)  also  predicted  a  second  Ha -emitting  disk  as¬ 
sociated  with  the  secondary  component.  Although  our  V 2  data  in 
the  Ha  channel  do  not  have  a  very  high  S/N,  we  do  not  detect  any 
deviations  from  the  elliptical  Gaussian  model  that  would  sug¬ 
gest  a  binary  signature  in  the  Ha  signal.  Therefore,  we  are  forced 
to  conclude  that  if  the  secondary  star  does  possess  a  Ha  -emitting 
disk  it  does  not  contribute  significantly  to  the  net  Ha  emission 
line.  This  is  also  consistent  with  the  secondary’s  disk  radius  of  6- 
8  R  .  estimated  by  Bozic  et  al.  (1995)  based  on  the  peak  sepa¬ 
ration  of  the  He  11  emission  line.  Assuming  that  the  Ha  emission 
is  proportional  to  the  area  of  the  disk,  as  demonstrated  by  Tycner 
et  al.  (2005),  we  conclude  that  the  Ha  emission  from  the  disk  of 
the  secondary  contributes  less  than  ~1%  to  the  total  Ha  flux. 

The  direct  detection  of  the  secondary  with  the  NPOI  is  un¬ 
likely  as  Gies  et  al.  (1998)  report  a  brightness  ratio  of  0.15  (or 
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A  m  of  2)  at  164.7  nm.  The  NPOI  operates  over  a  bandwidth  of 
550-850  nm,  and  at  these  wavelengths  the  magnitude  difference 
between  the  two  components  is  likely  to  be  much  larger  than  2.  A 
search  for  a  binary  signature  (a  sinusoidal  variation)  across  the 
continuum  channels  in  both  7  Cas  and  (j)  Per  yielded  null  results. 
On  the  other  hand,  a  binary  signature  with  a  small  V2  amplitude 
(less  than  ~3%)  could  not  be  ruled  out.  In  fact,  it  might  be  pos¬ 
sible  to  detect  the  signature  of  binarity  in  both  stars  at  longer 
wavelengths,  since  7  Cas  may  be  associated  with  a  cool  compan¬ 
ion  and  (l)  Per  might  contain  larger  contributions  from  the  circum- 
stellar  disks  of  both  components  due  to  free-free  and  free-bound 
emission. 

6.  CONCLUSIONS 

We  have  successfully  demonstrated  the  use  of  a  narrowband 
filter  in  the  NPOI  to  increase  the  contrast  between  the  Ha- 
emitting  material  and  the  central  star.  Our  observations  of  two 
Be  stars,  7  Cas  and  0  Per,  and  our  subsequent  analysis  have 
yielded  the  following  results: 

1 .  We  have  demonstrated  that  the  unifonn  disk  or  ringlike  mod¬ 
els  are  inconsistent  with  the  observations  of  the  circumstellar  re¬ 
gion  of  7  Cas,  whereas  a  Gaussian  model  is  fully  consistent  with 
the  data.  However,  since  the  themial  structure  of  the  circumstel¬ 
lar  disks  of  Be  stars  can  be  quite  complex,  as  suggested  by  recent 
models  (Millar  &  Marlborough  1 998;  Carciofi  &  Bjorkman  2006), 
it  might  turn  out  that  a  more  complicated  brightness  distribution 
may  be  required  to  fully  describe  these  regions. 

2.  The  circumstellar  disk  of  7  Cas  appears  to  be  consistent 
with  the  orbital  parameters  published  in  the  literature.  However, 
higher  precision  binary  solutions  are  required  to  test  for  the 
possible  disk  truncation  by  the  secondary. 

3 .  Based  on  interferometric  and  spectroscopic  data  we  have 
shown  that  the  brightness  distribution  of  the  Her -emitting  cir¬ 
cumstellar  region  of  (j)  Per  is  best  represented  by  a  Gaussian 
distribution. 

4.  Our  analysis  supports  the  earlier  prediction  by  Waters 
(1986)  that  the  disk  of  (f)  Per  is  truncated  due  to  the  presence  of  an 
orbiting  companion.  However,  the  disk  size  we  report  is  different 
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than  the  value  determined  by  Waters  (1986).  This  discrepancy  is 
likely  due  to  different  wavelength  regimes  used  in  each  study 
and/or  due  to  the  low  inclination  angle  assumed  by  Waters  (1986) 
for  this  star,  whereas  we  obtain  i  <;  74°. 

A  natural  extension  of  the  analysis  presented  in  this  study  is 
the  direct  comparison  of  interferometric  data  with  theoretically 
predicted  interferometric  observables,  such  as  synthetic  squared 
visibilities.  Future  observations  and  modeling  are  planned  that 
combine  narrowband  interferometry  with  spectroscopy  for  sev¬ 
eral  other  Be  stars,  especially  those  that  possess  weak  Ha  emis¬ 
sion.  The  goal  of  this  work  is  to  place  greater  constraints  on  the 
spatial  extent  and  physical  properties  of  the  circumstellar  ma¬ 
terial.  These  constraints  will  allow  various  dynamical  models  to 
be  tested  with  greater  certainty. 
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